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We study the high density region of QCD within an effective model obtained in the frame of the 
hopping parameter expansion. The model still acknowledges the sign problem peculiar to non- 
zero chemical potential, but it permits the development of refined algorithms which ensure a good 
overlap of the Monte Carlo ensemble with the true one. We review the main features of the model, 
including the most explicit form of the resumed expansion, and present calculations concerning 
the dependence of various observables on the chemical potential and on the temperature, in par- 
ticular of the charge density and the diquark susceptibility, which may be used to characterize the 
various phases expected at high baryonic density. 
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1. Introduction 

The aim of this work is to understand the phase structure of high density, strongly interacting 
matter. In the spirit of the /i = quenched approximation one has considered a "non-zero density 
quenched approximation" for /J. > based on the double limit M —> oo^ n — > oo ? £ = exp (fi — InM) : 
fixed [[lj, ||]. This implements a static, charged background, which influences the gluonic dynamics 
[Q, §1. The present model [Q] represents a systematic extension of the above considerations: the 
gluonic vacuum is enriched by the effects of dynamical quarks of large (but not infinite) mass, 
bringing a large net baryonic charge. In [^] and in the present paper we explore the phase structure 
of the model, as a first step in understanding the properties of such a background. 

This model can be derived as 1/M expansion of QCD at large jU around the unphysical massive 
quarks point. However, it is more realistic to understand it as an approximation, whose justification 
relies on the predominant role of the gluonic dynamics. We want to understand how this dynamics 
is influenced by the presence of charged matter. This would allow, among other things, to study the 
effect of dense, heavier background baryonic charges on light quarks and hadrons. 

The main ingredient of the model are Polyakov-type of loops, capturing the effect of heavy, 
limitedly moving quarks. The model still has a sign problem, but due to the factorization of the 
fermionic determinant, it allows for local, refined algorithms and large statistics. 

2. The model 

The QCD grand canonical partition function with Wilson fermions at /I > is [||] : 

2r(P,K,n) = J [DU]e- s ^ u »J2r F (K,n,{U}), 3f F {K,n,{U}) = DetW(K,fi,{U}), 

3 

W ff = 8 ff [\ - K f £ (T+iUiTi + r^TfU*) - Kf (e^r^Ti + e-^rUW)], 



r ±u = l±y u , k=— - = — -, (2.1) 

M mi 2(M + 3 + coshjU) 2(M + 4)' v 

with M the "bare mass", Mq: bare mass at jU = 0, /: flavour, U : links and T: lattice translations. 
The hopping parameter expansion is an expansion in closed loops of DetW = fT/DetW/, 



DetW/ =exp 



oo 4 

nnn^ 1 - (*ma^ 

l=Y{%}a=\ 



l=\ms=\ 

(2.2) 

g f % = (ee ±N ^f) r if % : "Polyakov r-path" , = 1 else ^% = \\T U = \\U l (2.3) 

where % are distinguishable, non-exactly-self-repeating closed paths of length /, s is the number 
of coverings of %, and a "Polyakov r—path" winds r times over the lattice in the ±4 direction 
with periodic(antiperiodic) b.c. e = +1(— 1) (p.b.c. in spatial directions). The factors are the 
eigenvalues of r 3^g l and the expansion can be expressed as a product of colour determinants p]. 
Notice that DetW is a polynomial in K of the order An c rifNy (colours, flavours, lattice volume), 
which means that the expansions on the right hand side can be truncated at this order. 
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Figure 1: Loops in the determinant (left); quark (middle) and di-quark propagators (right). 



The "quenched limit at /f > 0" is [g]: 

k^0,v^°°, Ke^ = i;: fixed, C = (l£f x 
^ 0] (C,{U})=exp 

L W*=i 

and the next order corrections lead to M] 



££ (£Cf Tr<; 



nDetc(l - £C^)\ 

{3} 



2f®(K,H,{U}) = exp^ -2£ £ ^Tr 



£C)' S 



{2} ,= 1 

(C, {t/}) n Det c (l - (SCY K 2 0%^ 



x,r,q,i,t,t' 



r,q,i,t,t' 

2 



(2.4) 
(2.5) 



(2.6) 



See Fig. |]. For easy bookkeeping we use the incomplete temporal gauge: 

U nA = 1 , except U { ^ n4=Nt)A = Vj, (2.7) 
&3 = V 3 , 9% f j = (V s y-W mi (V- + i)«U^ tl)i , r > q > 0, ±i= 1,2,3, 1 < t < t' < (2.8) 

(t < t' for q = 0). Eqs. ( p4[^8| ) define the present model. Notice that for U £ SU(3) we have: 

Det c (l - eCU) = 1- sCTrU + ^Trll^ -EC 3 (2.9) 

We measure spatial and temporal plaquettes, Polyakov loops and baryon charge densities 



r3"3AT3"' f/ 5 --^^y g 



(2.10) 



as well as diquark propagators (in maximal temporal gauge) and the corresponding susceptibility, 
which could give information on the large ji. physics [Jsj] (see [||] for explicite formulae, see Fig. [I]). 



3. Algorithm and Simulation 



We use the Wilson action and Wilson fermions within a reweighting procedure. The updating 
is done with a local Boltzmann factor (which only leads to a redefinition of the "rest plaquette"): 



J-J jTrPlaq "Q exp J 2 CMt C 



Plaq 



(3.1) 
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We here employ the Cabbibo-Marinari heat-bath procedure mixed with over-relaxation. This up- 
dating already takes into account part of the pL > effects and the generated ensemble can thus 
have a better overlap with the true one than an updating at /i = 0. Notice that such a factor may 
also improve convergence of full QCD simulations at jj. > 0. The weight (global, vectorizable) is: 



ri ex p 



-2C^eTr ( 



x,i,t,t' 



(3.2) 



One also can use an improved partition between the updating factor and the weight, to be taken 
care of by a supplementary Metropolis check. Anisotropy can be straightforwardly introduced. 

The simulations are mainly done on lattices 6 4 and 8 4 (in the incomplete temporal gauge (2.7)). 
We shall present results for tif = 1 and tif = 3 degenerate flavours (any mixture of flavours can be 
implemented). The fc dependence has been analyzed in [ffl]. Here we set K = 0.12, which drives the 
1 /M 2 effects in the baryonic density to about 50%. Our problem setting is primordially to explore 
the phase structure of the model at large jx, small T and we accordingly vary j8 and /x. 



4. Results 

Some analytic results can be obtained in strong coupling. To lowest order in j8 we have: 
1 



3 lr W^ *' 3 C 1+4C 3 



1+4C 3 + C 6 
C 3 

+c 6 



2N T „ 2 1 + 3C 2 + 2C 3 + 2C 5 + C 6 

1+ — " r ttw 

2N T o -2 1 + 3C2 + 2C3 + 2C5 + C6 



1+C 3 



(4.1) 
(4.2) 



for rif = 1. Fig. g shows increasing deviation of the full results from strong coupling with increas- 
ing j8 and indication of qualitative change of behaviour at large /I and above /3 =5.5. 
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Figure 2: Polyakov loops P (triangles) and (diamonds) vs C, 4 4 (left, j3 = 3, full and j3 = 5, open 
symbols) and 6 4 (right, j8 =5.5, full symbols and /3 =5.6, open symbols), strong coupling (lines). 

On Fig. [| left plot, we show for tif = 3 the behaviour with j3 at fixed jit values, which in 
the physical T, pL plane would mean moving on lines T ~ ^rMphys- We see at both values of /i 
qualitative changes of behaviour suggesting transitions from low to high temperature phases at 
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values of /3 as indicated by the Polyakov loop susceptibility (curves on the plot). This appears 
corroborated by the behaviour shown on the right hand plot, where we vary }X at two fixed j3 values 
and see a delayed onset of the "transition" for the smaller /3. The diquark susceptibility (we only 
show the second order term) is rather flat in j3 and shows a signal only at jj. approaching 1 . 
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Figure 3: Baryonic density (triangles), Polyakov loops (circles) and diquark susceptibility (squares) vs /3 
(left, = 0.85 open, and ji = 0.90 full symbols) and ji (right, j3 = 5.5 full and f5 = 5.6 open symbols). 

The dependence on rif is illustrated on Fig. ||, both for the 6 4 and the 8 4 lattices. As one can 
see, the highest values of /i attainable in this method are about 5 ~ 6T, depending on n/. 
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Figure 4: Baryonic density (triangles) and Polyakov loops (circles) vs jX at j3 = 5.6, 6 4 (left) and 8 4 (right; 
here also the diquark susceptibility, squares): «/ = 1 open and «/ = 3 full symbols. 



5. Discussion 

Our calculations show strong effects at large fi, which even at moderate j3 depart considerably 
from strong coupling estimates and also indicate possible phase transitions. The results concerning 
the latter for the n/ = 3 case are summarized in Fig. g. On the left plot we indicate the lines we 
have followed in the simulation in the r,/i p h ys plane. On the right hand plot we show the /3,/i 
plane together with the points at which we see indication of transitions. Since we cannot yet fix 
a scale (but see below) we preferred to show these points in the bare parameter plane. While the 
indication for a transitions are clear, it is uncertain what happens at larger /J. and small T and how 
the transition lines run. This region needs therefore further study. 
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P i (T) 




Figure 5: Prospected phase structure (only indicative!) and observed changes of behaviour (transitions?). 



The algorithm works reasonably well over a wide range of parameters and for lattices up to 8 4 . 
We reach large densities and j for temperatures ~ ^ T c (ji. = 0) or less. It appears difficult, however, 
to go to larger lattices and larger /» with this algorithm and one should consider improving it. 

The model permits to vary /», K and T as independent parameters. It is therefore interesting to 
extend the study to cover this full variability. For higher orders in K the bookkeeping soon becomes 
unmanageable, one could however consider using statistical ensembles of large loops 

Concerning the significance of this analysis we can take two points of view: 

Firstly, we can consider this model for itself, as describing "quasi-static charges" interacting 
via gauge forces. One may then ask whether this dynamics may lead to a non-trivial phase structure. 

Secondly, we can consider this model as an evolved "quenched approximation" in the presence 
of charged matter. Then this study would give us information about the structure of the so modified 
gluonic vacuum of the SU(3) theory. It would then be natural to think of it as providing a heavy, 
dense, charged background for light quarks propagation and calculate light hadron spectra and 
other hadronic properties under such conditions. This could also help fixing a scale. 

Acknowledgments: The calculations have been done on the VPP5000 computer of the Uni- 
versity of Karlsruhe and on the PC Cluster of the Institute of Physics of the Parma University. 
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